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In this paper we use non-Gaussian hydrodynamics to study the magnetic response of a flux-line 
liquid in the mixed state of a type-II superconductor. Both the derivation of our model, which goes 
beyond conventional Gaussian flux liquid hydrodynamics, and its relationship to other approaches 
used in the literature are discussed. We focus on the response to a transverse tilting field which 
is controlled by the tilt modulus, C44, of the flux array. We show that interaction effects can 
enhance C44 even in infinitely thick clean materials. This enhancement can be interpreted as the 
appearance of a disentangled flux-liquid fraction. In contrast to earlier work, our theory incorporates 
the nonlocality of the intervortex interaction in the field direction. This nonlocality is crucial for 
obtaining a nonvanishing renormalization of the tilt modulus in the thermodynamic limit of thick 
samples. 

PACS: 74.60-w, 74.60Ec 



I. INTRODUCTION 

The static and dynamical properties of magnetic flux 
lattices in type-II superconductors have been the focus of 
much-|theoretical and experimental work over the last ten 
yearsM. Interest in this field was revived by the discovery 
of the high-T c materials, where thermal fluctuations melt 
the Abrikosov flux lattice at temperatures and fields well 
below the mean field transition at H C 2(T]au. The flux 
lattice melting is a first order transition in clean samplesu, 
with an associated jump in the bulk magnetization, and 
it has been observed experimentallylilj. In conventional 
low-temperature type-II superconductors, the region of 
the phase diagram where thermal fluctuations are impor- 
tant is extremely small and mean field theory provides a 
good description of the physics of the flux-line array. In 
the high-T c materials, in contrast, the melted flux liquid 
replaces the Abrikosov lattice over a large region of the 
phase diagram. Understanding the properties of the flux 
liquid is therefore crucial for controlling the magnetic re- 
sponse of these materials. 

The conventional Abrikosov flux lattice is character- 
ized by two broken symmetries. First, the translational 
symmetry is broken by the ordering of the magnetic flux 
lines in a triangular lattice in the plane perpendicular 
to the external field. Secondly, the gauge symmetry 
along the field is broken by the alignment of the vor- 
tices with the external field. A natural question then 
arises of whether these two symmetries are recovered si- 
multaneously upon melting, or rather they are recovered 
in succession at two different temperatures. The latter 
scenario would allow for the appearance of a disentan- 
gled flux liquid phase where translational symmetry is 
recovered, but the longitudinal gauge symmetry is still 
broken. At a second transition temperature the disen- 
tangled flux liquid would then be replaced by an entan- 
gled flux liquid where the longitudinal gauge symmetry is 
also recovered. Alternatively, if both symmetries are re- 
covered simultaneously, the Abrikosov lattice would melt 



directly into an entangled flux liquid. The precise natuxc 
of such an entangled liquid remains an open questionLl. 
The existence of a disentangled liquid phase, exhibiting 
longitudinal superconductivity - the ability to support 
currents flowing without dissipation in the direction par- 
allel to the flux lines - in clean samples has been proposed 
some time ago by Feigel'man and collaboratorsEJ. Effty 
simulations provided support for Feigel'man's ideasEjtfl, 
but more recent numerical work indicates that the two 
transitions observed in earlier woj'k,may have been the 
consequence of finite size effectsLLSEl Recent numerical 
results support the scenario that the Abrikosov lattice 
melts directly into an entangled liquid and no disentaii- 
gled liquid phase exists in infinitely thick samplesEJaEII. 
Open questions, however, remain concerning the role of 
various approximations used in the different numerical 
models, particularly the range of the intervortex interac- 
tion. 

A closely related property of the vortex array that pro- 
vides a direct measure of longitudinal vortex correlations 
is the tilt modulus, C44. It can be probed by measur- 
ing the response of the flux array to a small additional 
magnetic field 6H±, applied perpendicular to the exter- 
nal field zHq responsible for the onset of the vortex state. 
Such a transverse field tilts the lines away from the direc- 
tion of alignment with Hq. Correlated disorder induced, 
for instance, by aligned damage tracks in the material 
can drive I/C44 to zero, yielding a transverse Meissner 
effect, which has-been proposed as the signature of the 
Bose glass phaseEaEi The role of correlated disorder in 
enhancing C44 in the liquid phase has also been observed 
experimentally in materials with a single family of twin 
planes by using the dc flux transformer configuratioro. 
These materials contain practically no small-scale disor- 
der, so that the macroscopic flux liquid regions in the 
channels between twin planes are very clean. The ex- 
periments suggest that the enhancement of C44, inter- 
preted as the onset of disentangled liquid phase, be a 
finite-size effect, that decreases with increasing sample 
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thicknesses. In thick samples the experiments indicate 
that the vortex-lattice melting and the loss of longitudi- 
nal superconductivity coincide in clean materials. Even 
though a true Meissner effect with vanishing 1 /C44 is not 
expected in infinitely thick, clean samples, it is clear that 
interactions can enhance the tilt modulus of clean flux 
liquids and suppress the transverse response of the su- 
perconductor. 

In this paper we employ hydrodynamics to evaluate 
the renormalization of the tilt modulus of a clean flux 
liquid due to interactions. Our starting point is a long- 
wavelength hydrodynamic free energy that includes non- 
Gaussian couplings in the hydrodynamic fields. It there- 
fore goes beyond the Gaussian flux-line liquid hydr. 
namic free energy discussed before in the literatureE£ll 
We show that such a non-Gaussian hydrodynamic free 
energy can either be written down phenomenologically 
or it can be derived by using the mapping of the clas- 
sical statistical mechanics of vortex lines with nonlocal 
interactions onto the quantum statistical mechanics of 
two-dimensional charged bosons, introduced some time 
ago by Feigel'man and collaboratorsEij. Our central re- 
sult is the expression for the renormalized wave vector- 
dependent tilt modulus given in Eq. 1.7 below. This 
is a perturbative. result that extends earlier results by 
other authorsE3 : E2l in two important ways. First, it incor- 
porates both the finite range and the nonlocality of the 
intervortex interaction in the field direction. This nonlo- 
cality plays a crucial role in controlling the tilt response. 
It is only when the nonlocality is properly accounted for 
that a finite renormalization of C44 is obtained in a clean 
flux-line liquids of infinite thickness. In addition, our 
formalism allows us to evaluate the full wave vector de- 
pendence of the renormalized tilt modulus - a result that 
was not discussed before in the literature. 

Before discussing our result in more detail, it is use- 
ful to make contact with already existing work. The tilt 
modulus of the Abrikosov lattice is easily calculated from 
the Ginzburg-Landau free energy for a superconductor in 
a field. It is dispersive both in the longitudinal and in the 
in-plane directions due to the non-local character of the 
intervortex interaction and it has a rather complicated 
expression, particularly for layered material. It naturally 
separates in the sum of two contributions, 



C44(<5LL, q z ) = c v i4 {q z ) + c^(q±, q z ), 



(1.1) 



with q± and q z wave vectors perpendicular and parallel 
to the external field, re spe ctively. The first term on the 
right hand side of Eq. (1.1) is the single vortex contribu- 
tion, arising from the self-energy part of the tilt energy. 
Neglecting its weak logarithmic dependence on q z , it is 
given byEHd 



S44 ~ n o£l, 



(1.2) 



where n$ = Bq z /(/)q is the average areal density of vor- 
tices, with Bq z the mean induction along the external 
field direction and </>o — hc/2e the flux quantum, and ei 



is the single- vort ex tilt energy defined below. The second 
term in Eq. (1.1) is the compressional contribution from 
intervortex interactions. It is stratLgly dispersive and in 
layered materials it is given by 
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where = Aj_/(1 — H/H^) 1 ^ 2 is the effective pene- 
tration length in the ab plane (the field is applied along 
the c axis) and p is the anisotropy ratio. It is important 
to stress that the long wavelength tilt modulus, 



c 44 = c 44 (q ± = 0,q z = 0) 
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is generally dominated by the large compressional con- 
/47r). The second term inside the brackets 



tribution (Bq z 



in Eq. (1.4), arising from the single- vortex contribution, 
is important only at very low vortex densities. 

The tilt modulus of a flux-line liquid cannot be eval- 
uated directly. It is, however, expected that the bare 
flux-liquid tilt modulus, denoted here by C4 4 (qj_, q z ), does 
not differj-Cjonsiderably from that of the lattice given in 
Eq. ([Ll])E3. In fact, a direct coarse-graining of the mi- 
croscopic intervortex interaction yields a Gaussian long- 
wavelength free energy of an entangled flux-line liquid 
with a tilt modulus given precisely by Eq. (1.1) aboveE3. 



Interactions responsible for nonlinearities in the long- 
wavelength free energy will, however, renormalize C44. 

The renormalization of C44 in flux-line liquids has been 
studied before by employing the analogy between the di- 
rected vortex lines induced in a three dimensional super- 
conductor by the external field zHq and thejmaginary- 
time world lines of two dimensional bosonsMlirQ. The 
most severe approximation made in-,the form of this bo- 
son mapping introduced by Nelsonuu, is that the pair- 
wise interaction between flux lines is approximated as 
local in the field direction (z), i.e., only the interaction 
between vortex segments at equal height z is considered. 
This corresponds to an instantaneous pairwise interac- 
tion between the bosons. One of the consequences of this 
approximation is that it completely neglects the com- 
pressional part of the tilt modulus. Hence in this model 
C44 is given by the single vortex part, which is inversely 
proportional to the boson superfluid density, n s , 
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The superfluid phase of bosons (n s = no) corresponds to 
an entan gled liquid of magnetic flux lines with C44 given 
by Eq. (1.2). A finite normal-fluid fraction of bosons 
of density n n = no — n s corresponds to a disentangled 
fraction of flux liquid and enhances the tilt modulus. A 
normal-fluid phase of bosons with n s = corresponds 
to a disentangled flux liquid with infinite tilt modulus 
and transverse Meissner effect. Tauber and Nelson (TN) 
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recently employed this boson mapping to evaluate the 
renormalization of c| 4 due to sample thickness, diffpft- 
ent boundary conditions and various type of disordered. 
They found that for finite sample thickness (correspond- 
ing to a nonzero boson temperature) there is a nonvan- 
ishing normal- fluid component which suppresses C44. On 
the other hand, the normal-fluid density always vanishes 
for infinitely thick samples (or vanishing boson tempera- 
ture), so that the flux liquid is always entangled in this 
limit. _ 

Feigel'man and coworkersEj incorporated the nonlocal- 
ity of the intervortex interaction in the field direction 
in the boson formalism. They showed that the statisti- 
cal mechanics of vortex lines with nonlocal interactions 
maps onto that of two-dimensional charged bosons. This 
nonlocal mapping incorporates the compressionaL part 
of the vortex tilt .modulus. Larkin and VinokureJ and 
later GeshkenbeinEJ used this nonlocal boson mapping 
to generalize the expression (fTj|) obtained by TN. These 
authors proposed that the long-wavelength renormalized 
tilt modulus can be written in terms of the superfluid 
density n s of two-dimensional bosons interacting with a 
gauge field as 
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The superfluid density was_evaluated perturbatively by 
Feigel'man and coworkersEj for the case where the re- 
pulsive interaction among the bosons is infinitely long- 
ranged, corresponding to a vortex liquid with A^ — > 00. 
These authors argued that in this limit a distinct dis- 
entangled flux liquid phase with diverging C44 exists in 
infinitely thick superconducting samples. 

The calculation of the interaction-renormalization of 
the flux liquid tilt modulus via hydrodynamics described 
here has the advantage that it naturally incorporates the 
nonlocality of the intervortex interaction and it allows 
us to easily treat the case of finite A. The non-Gaussian 
hydrodynamics used as the starting point contains bare 
elastic constants that are determined by the intervortex 
interact ion. In particular, the bare tilt modulus is given 
by Eq. (1.1). The corrections to C44 due to the nonlinear- 



ities are evaluated perturbatively. Our main result is an 
expression for the wave vector-dependent renormalized 
tilt modulus, given by 
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where n n (qj_, q z ) has the rather complicated integral ex- 
pression given in Eq. (6.7) below. The corrections to 
the tilt modulus incorporated in n n can be interpreted 
in terms of a disentangled fraction of the flux liquid - 
hence a "normal-fluid component" . When the nonlocal- 
ity of the intervortex interaction in the field direction is 



neglected, Eq. ( [L7| ) becomes identical to the result ob- 
tained by Tauber and Nelson (see Eq. (3.33) of Ref. p8| ). 
In this case the long-wavelength C44 is not renormalized 
in infinitely thick sa mple s. 

Our result, Eq. (1.7), is also simply related to the 
Larkin- Vinokur formula given in Eq. fly]). This is im- 
mediately seen by introducing a normal fluid fraction in 
Eq. (E6) as n n = hq — n Sl and then expanding for small 
values of the normal fluid fraction, n n /no << 1, to ob- 
tain 
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with C44 given by Eq. (jT|). This expression is formally 
identical to the long-wavelength (q± — 0, q z — 0) limit of 
our result. 

We find that interaction effects in a clean flux liq- 
uid do lead to a nonvanishing renormalization of the 
tilt modulus in the thermodynamic limit of thick sam- 
ples. This correction is present only if the nonlocality of 
the intervortex interaction is properly incorporated. The 
correction remains, however, small at all but very low 
(B < lTesla) fields. Our results are perturbative and 
cannot be used to infer quantitative conclusions about 
the existence of a true disentangled flux liquid phase. 
One of the main outcomes of our work is the development 
of a transparent hydrodynamic framework that can be 
used to study the role of the nonlocality of the intervor- 
tex interaction on the tilt response, both in clean materi- 
als and in the presence of disorder of various geometries. 
Note that in conventional, Gaussian hydrodynamics the 
effect of disorder on C44 cannot be detected. 

In section II we discuss the general form of the London 
free energy used as the starting point to study the mag- 
netic properties of superconductors in the mixed state. 
The various response functions of interest are also defined 
there. After discussing the response to a tilt field in sec- 
tion III, we review and contrast in sections IV and V, re- 
spectively, the results obtained by conventional Gaussian 
hydrodynamics and by the local boson mapping. After 
showing how hydrodynamics can be derived from the bo- 
son model in section VI, we introduce our non-Gaussian 
hydrodynamic model and discuss its relationship to pre- 
vious work. Our results are discussed in section VII. Fi- 
nally, a rigorous derivation of the nonlocal, non-Gaussian 
hydrodynamics from the charged boson analogy is dis- 
played in Appendix A, and the perturbative evaluation of 
the renormalization of C44 from interactions is displayed 
in Appendix B. 



II. MAGNETIC RESPONSE OF THE VORTEX 
ARRAY 

High-T c superconductors are uniaxial, strongly type-II 
materials with very large values of the Ginzburg-Landau 
parameter k = A/£. For applied fields H c i « H « 
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H C 2, their mixed state can be described in the London 
limit with a Ginzburg-Landau Hamiltonian given by 



H[0,A] 



4^A 2 V2vr M 



-in 



(V x A) 



(2.1) 



Here the z direction has been chosen along the anisotropy 
(c) axis of the superconductor. Greek indices (J,, v, ... run 
over all Cartesian components (/i = x, y, z) and summa- 
tion is intended in Eq. ( [2.l| ). Latin indices i,j,k, ... run 

only over x and y. The integral f t ... = f^dzf dr±... is 
over the volume f2 = LA of the superconductor, with L 
the thickness in the direction of the c axis and A the area 
in the ab plane. Also, A M = A M /(1 — H/H^) 1 ^ 2 , where 
A x = X y = Aj_ are the penetration depths from supercur- 
rents in the ab plane, while X z = p\± is the penetration 
depth from supercurrents along the c axis, with p the 
anisotropy ratio arising from an effective mass tensor for 
the superconducting electrons (p = [m z jm\}f l l' 1 ). Fi- 
nally, A is the total vector potential, with B = V x A 
the internal field in the material, and (/>q = hc/2e is the 
flux quantum. The corresponding Gibbs free energy func- 
tional is 



ones. We will specifically consider situations where the 
magnetic field responsible for the onset of the vortex state 
is applied along the z direction. Vortex line configura- 
tions are then conveniently characterized by a set of N 
single- valued functions r n (z), which specify the position 
of the n-th vortex line in the xy plane as it wanders along 
the z axis. The three-dimensional position of each flux 
line is parametrized as R n (z) = [r„(z), z] and the vortex 
density vector can be written as 



JV 



t(r) = £ rfRn(f) 5 ( 2)(r| 



dz 



r„(z)), 



(2.4) 



where r = (r±,z). The vortex density vector can be 
written as T(r) = (t,n), where t is a two-dimensional 
vector describing the local tilt of flux lines away from the 
direction of the external field and n is the areal density 
of vortices, 



n(r) = ^ 2 >(r x -r„(*)), 

71=1 

t(r) = E%^(r,-r„(z)). 



n=l 



(2.5) 
(2.6) 



0[0,H] =H[0,A] 
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(2.2) 



where H = V x A oxt is the applied external field. 

The London free energy functional can be rewritten in 
terms of interacting vortex lines by introducing a "vortex 
line density" vector defined as 



The vortex density vector is also directly related to the 
superfluid velocity of the electrons in the superconductor, 
v s = (fo/2n)V6- A, by 



6 T - B = V x 



(2.7) 



T(r) = x (V6 



(2.3) 



Here and below a hat (") is used, when needed, to dis- 
tinguish microscopic fluctuating quantities from average 



The Cartesian components of the local supercurrent are 
j s = (c/47rA 2 1 )w^ (no summation over [i intended here). 
After some manipulations (see, for instance, Ref. [l7] for 
the details) and neglecting spin wave fluctuations, one 
obtains 



S[T,H] = ^E{ [W(q) - 4(q)]tV(q)[<A)t(-q) - M^)} + |B(q)| 2 - 2H(q) ■ B(-q) j, 



with 



X±q 2 



A^i + Aig 2 2 



(2.8) 



(2.9) 



Here, q = (q_i_,<7 z ) and P?(q±) — Sij — q±iq±j is the two-dimensional transverse projection operator, with qj_ 



C l±/<1±- The corresponding longitudinal projection operator is P£(qjJ 



In this paper we will only consider magnetic field fluctuations due to fluctuations in the vortices' degrees of freedom. 
This London part of the field fluctuations is obtained by minimizing the Ginzburg-Landau free energy ( |2.8| ) for fixed 
vortex configurations T(q) and it is given by 



B(q)=B^(q)+B"(q) 
where B v (q) is the part of the internal field due to the vortices, 
Sj(q) = (1 + U^r'tMcOtof^q) 



(2.10) 
(2.11) 
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1 + X\q 2 
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1 + AiQ 2 



Wq), 
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and B (q) is the Meissner response of the material to a spatially inhomogeneous external field, 



Bf(q) = (l + U(q));X(q) 
1 



(2.12) 



1 + \ 2 ± q 2 



Ai<zV + ( 1 Al) g 2 Pg(qx) 



1 + Aig2 + A 2 g 



^(q)- 



In addition to the contributions given in Eq. (2.10), there are field fluctuations re prese nting therm al deviati ons from 



the solution of the London equation, which are neglected here. By inserting Eqs. ( 2.11 ) and ( 2.12 ) into Eq. (2 
obtain the vortex free energy functional expressed entirely in terms of vortex degrees of freedom, 



, we 



S[T,H] 
where 



^^{^(q)V( q ):t(-q) - ^(q^q^-q) - ^^(q)(l + U(q)) M >,(-q)|, (2.13) 



V(q) = ^o(i + u( q )) M V CT -( c i) 

V 



(2.14) 



i + \W 



S 5 5 ( A '~ A i)gj p T ( Ql ) 



are the Fourier components of the anisotropic intervortex interaction, with Vq = c/^/Att. One important property of 
the intervortex interaction is its nonlocality. In particular, the nonlocality in the z direction, reflecting that flux-line 
elements at different z heights repel each other via a Yukawa-like potential, will play a very important role in the 
discussion below. 



The Gibbs free energy of the vortex system is given by 
G(n,T) = -k B T1n.Z(IL,T), (2.15) 

where 



homogeneous field H = zHq and the local field is to be 



Z(H,T) = J Vf{r)e- g/kBT 



(2.16) 



is the canonical partition function. The prime over the 
integral sign indicates that the integration must be per- 
formed with the constraint V • B = 0. The average local 
field in the superconductor is then given by 



understood as the field given by Eq. (2.11). 

The focus of this paper is on the response of the vor- 
tex array created by the external field zHo to a small 
additional spatially inhomogeneous external field SH(r). 
The Gibbs free energy functional in the presence of this 
perturbation can be written as 

G(T, zH + SH) = g (T, H ) + 8G(T, 5H), (2.19) 



where Go is given by Eq. (2. IS) and the perturbation is 



B(r) =< B(r) >= -4tt 



SG 



(2.17) 



Sg(T,SB.) = -i- / B v SH 

47T 



(2.20) 
(2.21) 



where the brackets denote a statistical average with 
Boltzmann weight ~ exp[— Q/ksT]. 

For a spatially homogeneous external field applied 
along the 2- direction, H(r) = zH , we obtain the fa- 
miliar formlil, 



-NL 



H (j>o 

47T 



(2.18) 



^^(^^(q^C-q). 



2Q, 



For a uniform applied field H = zHq, the Meissner part 
of the transverse local field given in Eq. ( 2.1 2| ) vanishes. 
The local field in the superconductor is entirely due to 
the vortices and it given by Eq. (2.11). From here on 
we will always refer to the vortex system created by the 



The local field B v in Eq. Q2.20] ) is the field in the ab- 
sence of the perturbation SH and i s related to the vortex 
degrees of freedom via Eq. ( 2.11 ). It does not include 
the Meissner response to the perturbation 511. The su- 
percurrent is defined as j s = (c/47r)V x B y . 

Below we will use (...)o to denote a statistical average 
over the unperturbed ensemble described by Go, while 
(■■■)h wih denote the average over the pert urbed ensem- 
ble, with free energy given by Eq. (2.19). The mean 
local field TS H in the material in the presence of the per- 
turbation (5H can be written as the sum of vortex and 
Meissner parts as 



B H ( q ) = (B v (<i)) H +<5B M (q), 



(2.22) 



where (5B M (q) is the Meiss ner response to the pertur- 
bation, given by Eq. (Eu|) with H(q) = <5H(q). To 
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linear order in the perturbing field, the vortex contribu- 
tion can be expressed in terms of correlation functions in 
the unperturbed ensemble as, 



{B^(q)) H = (S^(q)) 



A<£J(q)4Vq)>^(q) 



Finally, in order to make contact with the literature, 
it is useful to write the perturbing field in terms of a vec- 
tor potential, <5H = V x <5A ext . The linear response to 
the vector potential (5A cxt is then characterized by the 
helicity tensor T^, which relates the induced current to 
<5A ex \ 



(2.23) 



where (...) c is the connected part of the correlator, i.e., 
(AB) C = (AB) - (A)(B). Finally, the corresponding lin- 
ear response function defines the magnetic susceptibility 
Xij (q) °f the material according to 

<(q) - (<(q))o = [47rx^(q) + S^SH^q). (2.24) 

The components of the susceptibility tensor can also be 
expressed in terms of vortex density correlations, 



47rx ([1I/ (q) = ~f- 



Vq k B TV < 

where T^ u (q) is the correlation function of the vortex 
density vector, 



(2.26) 



The density-vector correlation function can be expressed 
in terms of derivatives of the partition function of the 
perturbed system as 

(T^q^q')^ = {fokBTfiV- 1 )^- 1 )^ (2-27) 
-S 2 \nZ(H i + 5U,T)- 



5H K (q)6H x (q>) 



<5H=0 



where (V^ 1 )^ are the components of the inverse of the 
interaction tensor (2.14). 

The tensor is block diagonal, with T M „ = (T^, T zz ). 
The component T zz is the density-density correlation 
function or structure function of the vortices, 



T zz (q) = S(q) = (Sn(q)Sn(-q)) 0) 



(2.28) 



where dn(q) = n(q)— riofW q o describes the fluctuation of 
the local density field from its mean value uq = Bo z /(j}Q, 
with Bq z « Hq the equilibrium value of the z compo- 
nent of the internal field. The in-plane part Tij is the 
tilt-tilt autocorrelator and it is the central quantity of 
interest here. It can be written in terms of transverse 
and longitudinal components as 

TiM = r z (q)J$(qj.) + T r (q)if (q x ). (2.29) 

The transverse part of the tilt autocorrelator determines 
the tilt modulus of the vortex array. The wave-vector- 
dependent tilt modulus is defined by 



T r (q) = 



Cii{q±,q z ) 



(2.30) 



j"(q) = -cT^(q)6AT(q), 



(2.31) 



where j H is the total screening current induced in the ma- 
terial by the perturbing vector potential, comprising of 
both the vortex and Meissner contributions. The helicity 
tensor can be immediately related to the components of 
the susceptibility tensor, 



(2.32) 



Using Eq. (2.25), it can also be expressed in term of the 
correlations of the vortex density tensor. 



2V(q)K A (-q)r ffA (q), (2.25) 



III. TILTING FIELD 

In the remainder of this paper we focus on the response 
of the vortex array to a spatially inhomogeneous field 
<5H^(q) applied normal to the direction of Ho and that 
tilts the flux lines away from the z direction. As dis- 
cussed by Chen and TeiteO, we distinguish two types 
of perturbations. The first is a tilt perturbation, corre- 
sponding to a tilting field which is spatially homogeneous 
in the xy plane and may be modulated in the z direction. 
The long wavelength response to this tilt perturbation is 
determined by the long wavelength tilt modulus, C44, de- 
fined as 



n\k B T 
C44 



= lim lim T T (q±,q z ). 

9z^0 q±^Q 



(3.1) 



The order of the limits (q± — > first, followed by q z — » 0) 
is important here and reflects the physical situation of the 
relevant experiment. The vanishing of the long wave- 
length tilt modulus signals the onset of a transverse 
Meissner effect, where the perturbing field is c omp letely 
expelled from the material (as seen from Eq. ( |2.13| ), the 
corresponding static susceptibility equals — 1/4tt). This 
occurs, for instance, in vortex arrays pinned by columnar 
defects. 

The second physical experiment of interest here is the 
response to a tilting field (qj_ ) which is spatially ho- 
mogeneous in the z direction (i.e., independent of q z ) 
and generates a shear perturbation of the vortex array. 
Such a field can be obtained from a vector potential 
<5A ext = zcL4° xt (rjJ, which induces screening currents 
along the z direction. In the literature the response of 
the superconductor to such a shear perturbation is often 
characterized by the corresponding compon ent o f the he- 
licity modulus (T 2Z (qj_)) defined in Eq. ( 2.31 ), which 
in turn is related to the transverse part of the tilt-tilt 
correlator by 



G 



4ir 1 
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1 - 



V T T (q ± ,q z =0) 



k B T 1 + q \\\ 



(3.2) 



where the first term arises from the Meissner part of the 
response. The long wavelength limit of the helicity mod- 
ulus is 



„2 r 



lim T zz (q±) = — 

q±^0 47T 



1 



V 



k B Tg^o 



lim T T (qj_,q z = 0) 



(3.3) 



The vanishing of lim q± ^ a Tx(q±,q z = 0) yields 
lim^^o AttT zz (q±) / q\ = 1, which corresponds to a per- 
fect Meissner response in the z direction and signals lon- 
gitudinal superconductivity. 

We emphasize, however, that both the perturbations 
just described simply probe the magnetic response of the 
superconductor, which is the true equilibrium test of su- 
perconductivity. In fact the relevant response function in 
each case (tilt or helicity modulus) is simply the trans- 
verse part of the susceptibility tensor, 

Xr(q) = if (qx)Xy(q)- (3-4) 
The long wavelength tilt modulus is given by 



n 2 V 
C44 



1 + 4tt lim . xr(q± = 0,q z ), 



(3.5) 



and the component of the helicity modulus that controls 
longitudinal superconductivity is 

lim T zz (q±)=- lim q± 2 xr(q±, q z = 0). (3.6) 

q±^0 q±^0 

In a flux-line lattice the transverse part of the tilt-tilt 
correlator is non-analytic at small wave-vectors and the 
different order of limits of the two perturbations dis- 
cussed above is important. This is because the vortex 
array has a nonzero long wavelength shear modulus, Cqq. 
As a result, the flux lattice exhibits longitudinal super- 
conductivity, with \im q± ^Q Tx(qj_, q z = 0) = 0, and 



limyT tfa (9i 

q±->o 



■ qz 



0) 



1 

"4tt' 



(3.7) 



but no transverse Meissner effect, as lim gz _>o 7Y(q^ = 
0, 9z )^0and 



lim XT ttice (gj 



0,g z 



1 

47T 
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(3.8) 



In a flux- line liquid, in contrast, we find that the order 
of limits is not important and the flux array in general ex- 
hibits neither longitudinal superconductivity, nor perfect 
Meissner effect, as 

lim Xr qU1C W = °) Qz) = lim Xr" 1 (qx , Qz = 0) (3.9) 



1 

47T 



'-44 



where cf~ 4 is the flux liquid tilt modulus, renormalized by 
interaction effects. We will see below, however, that in- 
teractions can yield a strong upward renormalization of 
C44 even in clean flux liquids. 



IV. GAUSSIAN HYDRODYNAMICS 



A useful framework for discussing the long wavelength 
properties of flux-line liquids that naturally incorporates 
all nonlocalities of the intervortex interaction is hydrody- 
namics, where vortex fluctuations are described in terms 
of a few coarse-grained fields. By long wavelengths, we 
mean wavelengths large compared to the spacing between 
CuC>2 planes in the z direction, and large compared to 
the intervortex spacing in the ab plane normal to z. 

The coarse-grained hydrodynamic fields for a flux-line 
liquid are the fluctuating areal density, 



JY 



H (r) = Y / S { B 2 l(r ± -r n (z)), 



71=1 



(4.1) 



and a tilt field, 



N , 

F ( r )-E^z(^-r„(z)). 



(4.2) 



Here S BZ (r±) is a smeared-out two-dimensional 5- 
function with a finite spatial extent of the order of the 
inverse of the Brillouin zone boundary ksz — y/^wno- It 
is defined as 



E • 

q±<k B z 



(4.3) 



We stress that these hydrodynami c fie lds di ffer from the 
microscopic fields defined in Eq. ( |2.5D and ( |2.6| ) as they 
are coarse-grained quantities obtained by averaging out 
the more microscopic and rapidly varying degrees of free- 
dom. 

A Gaussian hydrodynamic free energy containing 
terms quadratic in the deviations of the fields from their 
equilibrium values can be obtained by coarse-graining the 
microscopic energy ofj-interacting vortices given in Eq. 
Q2.19D , with the result!-], 



Fa 



1 

An Jt Jt' 



B{r-r')Sh H (r)Sh H (r') (4.4) 
+K(r-r')i H (r) ■ i H (r')}, 



where Sn H (r) = h H (r) — uq and B(r) and K(r) are non- 
local liquid elastic constants. The density and tilt fields 
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are not independent quantities, but are related by a "con- 
tinuity" equation expressing the constraint that vortex 
lines cannot start or stop inside the sample, 



and 



d z 8n 



H 



V ± • t 



H 



0. 



(4.5) 



The Gaussian hydrodynamic free energy is rewritten in 
a more familiar form by passing to Fourier space, 



F G = 



1 



2nf.fl ^ 
u q 



(4.6) 



where cj 1 (q) and c\^(c\) are the bare compressional and 
tilt moduli of the flux liquid. The compressional modulus 
is given by 



c?r(q) = Bl - 



1 + q 2 \\p 2 



(4.7) 



The bare tilt modulus is found to be to a good approxi- 
matio n identi cal tp the flux lattice tilt modulus given in 



Eqs. (|L1 



4^ (l + g2 A 2J(l + g 2 A 2^2p2 A 2j 

It modul 
aitical tp 



In this Gaussian approximation, the probability of a 
fluctuation is proportional to exp(— Fc/k B T) and aver- 
ages must be car ried out subject to the continuity con- 
straint, Eq. ( |4.5| ). The correlation functions of the hy- 
drodynamic fields are then immediately calculated and 
are given by 



(<5n ff (-q)«5n ff (q)} ( 



(tf(-q)^(q)) ( 



n\k B Tq\ 



nlk B Tq ±i q z 



(4.8) 



(4.9) 



{tf (-q)tf (q)) G = I*(q)P£(<u) + T°( q )P^( qx ) 



with 



ir(q) 



n\k B T 



(4.10) 



(4.11) 



/{E|[^] 2 ^E^)-^)i) 

n— 1 m^n 



rg(q) 



n\k B Tq\ 



l (q)^+ C l(q ) g 2 



(4.12) 



The long wavelength tilt modulus is determined by the 
trans verse part of the tilt autocorrelator, according to 
Eq. (111). To this Gaussian order it is t hen identically 
given by its bare value, C44, given in Eq. (1.4). Gaussian 
hydrodynamics does not allow for any renormalization of 
the tilt modulus, even in the presence of disorder. This 
is because a disorder potential couples to the flux-line 
areal density that, within a Gaussian theory, is in turn 
decoupled from the transverse part of the tilt field. In 
particular, this naive hydrodynamic theory does not de- 
scribe the possibility of a disentangled flux-line liquid, 
with a tilt modulus enhanced by interaction or disorder. 
In other words, Gaussian hydrodynamics is by definition 
a theory of entangled flux- line liquids. 



V. 2D BOSON MODEL 

Considerable progress in understanding the properties 
of vortex-line arrays has been made by employing the 
formal analogy between the classical statistical mechan- 
ics of directed lines in three dimensions and the quantum 
statistical mechanics of two-dimensional bosons. The ad- 
vantage of this approach is that it can incorporate inter- 
action effects accounting for localization or disentangle- 
ment of the vortices. The drawback is that this model, 
at least in its simplest implementation employed by Nel- 
son and coworkersBoE3n3, neglects the nonlocality of the 
intervortex interaction. We will show below that the non- 
locality of the interaction in the field (z) direction plays 
a crucial role in controlling the tilt modulus. 

In this section we briefly review the local version c 
boson mapping employed by Nelson and coworkers 
and the results obtained recently for the tilt modulus by 
Tauber and NelsonEj using this model. 

Neglecting the nonlocality of the intervortex interac- 
tion, the free energy of interacting vortex lines in a field 
H = Hqz + 5Hj_ given in Eq. ( p. 19 ) is approximated as 



1 



g({r n },H) = NL(H ^-e 1 



(5.1) 



where ei = ti/p 2 , with ei = eqItik the effective line ( |5.l| ). Two crucial approximations have been made in 



tension and eo = (</>o/47rA_i_) a characteristic energy 
scale. The nonlocality relating fields and vortex vari- 
ables has been neglected also in the last term of Eq. 



rewriting the general intervortex energy given in (2. 
the form (5.1). Firs t, th e leading elastic term in the self- 



energy part of Eq. (2.8) has been linearized, according to 
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teraction among different ffux lines has been replaced by 
an interaction Vj_(rj_) acting locally in each constant-z 
plane, given by 



8tt 2 A 



(5.2) 



with Kq(x) a modified Bessel function. Of these approx- 
imations the latter is the most severe, since it amounts 
to neglecting the q z dependence of the elastic constants - 
an approximation that strongly affects the tilt modulus, 
as we will see below. Letting Q({r n (z)},H) = piNL + 
J%({r n (z)},H), with fi = H ^-e 1 = <f> { H o -#ci)/4tt 
a chemical potential, the grand canonical partition func- 
tion of the vortex liquid can be written as 



oo N « 

Z gr (H) = £ J] / Vr n (z)e-^" 

N=0 ' n=0 J 



)/k B T 



(5.3) 



The integral in Eq. (5.3) is over all vortex line configu- 
rations. It has the form of a quantum-mechanical parti- 
tion function in the path integral representation for the 
world lines of N particles of mass e%, moving through 
imaginary time z and interacting with the repulsive pair 
potential V±_ (r± ) . The vortex model with this simplified 
interaction can therefore be mapped into a model of 2D 
massive bosons with instantaneous pairwise interaction. 
The mapping results in the following correspondences: 



Z <-> T 

L ?l/3boson 

€\ «-> m 
ksT <-> h 



(5.4) 



H, 



4tt 
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where /?boson — l/^s^boson is the inverse temperature of 
the bosons. The precise mapping of the grand canoni- 
cal vortex line partition function onto the Feynman path 
integral in imaginary time r of a gas of two-dimensional 
bosons requires the introduction o f a second quantized 
Hamiltonian cocrcanonding to Eq. (5.1) and is described 
in the literaturea'E3c3. Some care must be taken in deal- 
ing with the tilting field SH± which introduces velocity- 
dependent terms into the fictitious boson Lagrangian. 
One important difference between the flux-line array and 
the boson system is in the boundary conditions in the fic- 
titious time variable z. The mapping of the free energy 
( |5.l|) of vortex lines onto the "action" of two-dimensional 
bosons is exact only when one imposes periodic bound- 
ary conditions for the flux lines in the z direction, i.e., 
r n (L) = r„(0). In contrast the natural boundary con- 
dition for flux line would be free boundary conditions, 



corresponding to (^jf_L_ L = {^f) z=0 = ^' ^ s snown 
by Tauber and NelsonEEl, the choice of the boundary con- 
ditions does affect the tilt modulus of a finite-thickness 
sample. We will not, however, discuss this here as we are 
ultimately interested in infinitely thick samples. 

To complete the mapping, the grand canonical parti- 
tion function (5T) is first rewritten in a coherent-state 
path integral representation as 

Z gr (H)= J Th/>(r ± ,z) J Pffn^je-*'*' 11 '/^. 

(5.5) 

The boson "action" in the imaginary-time path integral 
is 



ip*(k B Td z 



{k B Tf 
2ex 



VI V (5.6) 



J dv' ± V ± (r±-r' ± )\^(r ± ,z)\ 2 \ij(r' ± ,z)\' 



+ 



and h(r) = (</>o/47r)5Hj_(r). The complex fields ip and tp* 
correspond to boson annihilation and creation operators 
in the second quantized Hamiltonian. It is convenient to 
rewrite these fields in terms of an amplitude and a phase 



1>(t x ,z) = yjn(r ± ,z)e i0 ^>*\ 



(5.7) 



The magnitude n(rj_, z) of the field ip corresponds to the 
fluctuating local boson density. The phase field 9 deter- 
mines the boson momentum density, 



g(rj_, z) = k B TnV ±9 . 



(5.8) 



Upon inserting Eq. ( p.7| ) into Eq. (5.6), the action can 
be written in terms of density and phase variables as 



<^ikBTfid z 9 

(k B T) 2 (V ± n) 2 , {k B Tf 



(5.9) 
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n 



2ei 



n{V ± 9) 2 



inn -Vic/ n 



+ / dr'_ L V±(r_ L -r'_ L )n(r±,z)n(r'±,z) 



where we have dropped surface terms that vanish for pe- 
riodic boundary conditions. 

Th e tilt -tilt correlator Ty(q) can be calculated using 
Eq. fl2.27D , with the result, 







rp fs k B T 



(n(q)) h =c 



k B T 

vie! 



p.p' 



n(q-pMp)n(-q-p')0(p')>h= 



(5.10) 



where the b rack ets denote an average over the full non- 
linear action (5Jj), evaluated at h = 0. 

To proceed, a standard approximation is to consider 
only small fluctuations of the fields from their mean val- 
ues. Letting 



ra(rj_, z) = n + 5n(r±,z), 



(5.11) 



and retaining only terms quadratic in the fields in the 
action, the corresponding Gaussian action in zero tilting 
field is given by 



S G [4>,i/)*;0] = J ^ik B TShd z e 

, (k B T) 2 (V ± 6n) 2 , (k B T) 2 



(5.12) 
n (V ± #) 2 



8ei n 2e% 

V±(r± - T , ± )Sh(r_ L , z)Sh(r' ± ,z) 

To Gaussian order the tilt autocorrelator is given by 



no/.. , n $k B T r ( fn k B T\ 2 



q±iq±j(0(ci±)Q(-ci±))G 

(5.13) 



where (...)g denotes an average over the Gaussian action 
(5.12). The correlation functions of the fluctuating fields 
are easily calculated within the Gaussian approximation, 
with the result, 



(5n(-q)5h(q)) G = 



(0(-q)<fa(q)>G 



n k B Tq 2 /ei 
^e B {qf/{k B TY 



ql + e B { q fl(k B TY ' 

£ie B { q f/{n q 2 {k B Tf) 
q 2 + e B (q) 2 /(k B T) 2 ' 

where 

£B{q±) _ [ n k B Tq± 2 V±(q±) ,k B Tq ± 2 



(0(-q)0(q)) G = 



k B T 



+ (- 



2ei 



(5.14) 



(5.15) 



(5.16) 



1/2 



(5.17) 



corresponds to the Bogoliubov spectrum of the two- 
dimensional boson superfluid. The quartic term in the 
Bogoliubov spectrum arises from the |Vj_n| 2 "kinetic" 
term in the action. To this Gaussian order of approxi- 
mation the tilt modulus is dispersionless and simply the 
bare part of the single- vortex contribution to C44, given 
by 



(5.18) 



as given in Eq. (1.2). By compa ring the correlation 
functions given in Eqs. (5.14-5.16 ) to tho se of the hy- 
drodynamic fields given in Eqs. (4.8- 4.10| ), we see that 
the results obtained by these two methods agree with 
each other provided we drop the term of 0{q\_) in the 
Bogoliubov spectrum (which is of higher order in the 
wave vector and therefore is consistently neglected in 
a long wavelength theory) and make the identifications 
C44(<7J_,<?*) = n ei and (^ 1 {q±,q g ) = n$V±(q±). The 
quantity that replaces the "Bogoliubov spectrum" in hy- 
drodynamics is a characteristic inverse length scale ^J 1 
that controls the decay of correlations along the z direc- 
tion, given by 





r*(?±)i 


) 


. k B T\ 



-1 1/2 



6s 1 (qj-,qz) = qx\ 



' Cn(g_L,g z ) 

C4 4 (<7_L,<7z) 



(5.19) 



Notice, however, that, in contrast to the boson spec- 
trum, the correlation length t; z depends on q z , not just 
on <7j_. This dependence arises from the nonlocality of 
the intervortex interaction in the field direction and will 
have important consequences on the renormalization of 
C44. Finally, we stress that the hydrodynamic tilt field 
does not simply map onto the momentum density of two- 
dimensional bosons, which in turn is related to the boson 
phase variable by Eq. 



(5.S) 



The boson momentum den- 
sity is to lowest order purely longitudinal while the tilt 
field always has a transverse part. 

Tauber and Nelson evaluated perturbatively the cor- 
rections to arising from terms beyond Gaussian in the 
free energycS. These corrections can be obtained by fac- 
torizing the fourth order correlator on the right hand side 
of Eq. (5.10) as_a product of Gaussian correlators using 
Wick's theorerrO. For the long wavelength tilt modulus, 
these authors obtained 



where 



n.„ = 



1 

„vR 



Lk B T 



1 



n ei 



n B 
n 



(27T) 



q± 



sinh 



Le B (q±) 
2k B T 



(5.20) 



(5.21) 



is the normal-fluid density of the two-dimensional boson 
liquid. The long-wavelength tilt modulus can also be 
written as 



1 

c 44 
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(5.22) 



10 



where nf = uq — i s the boson superfluid density. As 
easily seen from Eq. fT^l] ) and discussed in TNc3, the 
normal-fluid density is finite only for samples of finite 
thickness L, corresponding to a nonzero boson tempera- 
ture. In this case one obtains a renormalization of the tilt 
modulus due to finite-size effects. The sign of this cor- 
rection is sensitive to the choice of boundary conditions 
(the result for periodic boundary conditions is displayed 
here). The normal fluid density vanishes, however, for 
L — ► oo. The local boson model therefore predicts that 
the tilt modulus of an infinitely thick, clean supercon- 
ductor is unrenormalized and equals its bare value noei- 
In other words, the flux- line liquid is always entangled in 
the thermodynamic limit. 



VI. NON-GAUSSIAN HYDRODYNAMICS AND 
DISENTANGLED FLUX LIQUIDS 

Our goal in the remainder of this paper is to construct 
a non- Gaussian fully nonlocal hydrodynamic theory and 
use it to evaluate the renormalization of the tilt modulus. 
As a first step in this direction, in this section we derive 
a non-Gaussian hydrodynamic free energy from the local 



boson action given in Eq. (5.9). Of course such a hydro- 
dynamic theory neglects interactions that are nonlocal 
in z and will mainly be used as a guide for constructing 
a more general non-Gaussian nonlocal hydrodynamics in 
the next section. The non-Gaussian terms in the free 
energy renormalize the tilt modulus. When these cor- 
rections are evaluated perturbatively, the resulting C44 is 
identical to that obtained by by Tauber and Nelson us- 
ing the boson formalismcS. The main goal of this section 
is to emphasize the relationship between the boson for- 
malism and hydrodynamics and to stress that equivalent 
results can be obtained by either method. 

To derive the hydrodynamic free energy from the bo- 
son action, we-iemploy the method used by Kamien and 
collaboratorscHl for the formally analog problem of di- 
rected polymers in a nematic solvent. We begin by elimi- 
nating the term n(V^fl) 2 in Eq. ( 5.1C ) in favor of 
a new vector field P, via a Hubbard-Stratonovich trans- 
formation, with the result 



Z gr (H) 



where 



S'[P, n,0; h] = / \ ik B Thd z 6 + {k ^ )2 (V ^" )2 + %^nV ± • [k B TV - h] 
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(6.1) 



(6.2) 



~ [(k B TVf -h 2 }+ f VI (rx - r' x )n(r ± , z)n(r' ± , z) 
*£l Jr' ± 



If we integrate over P in Eq. (6.3), we return to the original nonlinear action. Instead we integrate over 9 which only 
appears linearly in the new action. This integration results in a ^-functional, yielding 



Z gr (H) = / VnVPexp 



/ ln(n(r)/n )le-^[ ft ^ h ]/ fcBT <5(a z Ti + V ± -^(fcBTP + h)), 



(6.3) 



with 



Hi 



^-h 2 + / V±(t± -r' ± )n(r ± ,z)h{r' ± ,z)}. 

Jr' 1 ' 



(6.4) 



In obtai ning E q. (6.4) we have discretized the nonlin- 
ear action (5.1C) in real space, according to 



/(r) -> vo^fi, 



(6.5) 



with wo an elementary volume, vq = ei/(^B-^ n o)' This 
is the volume of a box with base area equal to 1/no and 
height equal to the single-vortex entanglement length, 



2k B Tno 



(6.6) 



The term containing the logarithm of the fluctuating den- 
sity arises from the Jacobian of the functional integration 
over the full nonlinear action. It represents the nonlinear 
"ideal gas" part of the flux liquid free energy. 



Statistical averages have to be performed by integrat- 
ing over the fields n(r) and P(r) wi th t he constraint pro- 
vided by the ^-functional in Eq. ( |6.3| ). Com paris on of 
Eq. ( |3.3| ) to the hydrodynamic free energy (h6) of a 
flux-line liquid with the constraint (fh^) suggests a phys- 
ical interpretation for the auxiliary vector field P. The 
quantity n(k B TP + h)?i takes the place of the hydro- 
dynamic tilt field t H introduced in the previous section. 
The difference between the vector field P and the tilt 
field can be understood by-aioting that, as pointed out 
by Nelson and Le Doussa£3, the canonically conjugate 
momentum of the fictitious particle that corresponds to 
the n-th flux-line is p„ = i (ei ^ + h) . The vector field 
P can then be interpreted as a sort of "velocity" field, 
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while the tilt field i H represents the canonically conju- 
gate momentum density. The two differ in the presence 
of an applied transverse field h that contributes to the 
single- vortex "canonical momentum" . 

The relationship between the effective action S B and 
the hydrodynamic free energy of a tilted flux-line liquid 
is made more transparent by performing and additional 
change of variable that replaces the field P by a tilt field 



defined as 



t(r) = 4% B TP(r)+h(r)]. 



(6.7) 



The Jacobian of this transformation cancels the Jacobian 
of the Hubbard-Stratonovich transformation used earlier 
and we obtain, 



Z gr (H_) = / VhVie- SH ^' i '' h ^ kBT S(d z 



(6.8) 



with 



Sff[ri,t;h] = 



1 



2k B T 
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{k B Tf (V X A) 2 
4ei h 



h • t + / Vi(r_L - r' x )ri(r_L, z)h(r' x , z) 



(6.9) 



The effective action of a tilted flux-line liquid given in 
Equation ([3J^) becomes formally identical to the corre- 
sponding nonlinear hydrodynamic free energy, provided 
we make the identifications, 

n(r) <-> h H (r), (6.10) 
t(r) «-> F(r), 
no?! <-> c" 4 (q), 
n 2 V ± {q ± ) «-> c^q). 

The corresponding hydrodynamic free energy is nonlin- 
ear, but local in z, and it given by 

F £ [n H ,t H ;h} = k B TS H [n,i;h}. (6.11) 

The subscript indicates that only local interaction 
among the vortices has been retained in this hydrody- 
namic free energy. The free energy F l contains the term 
quadratic in the density gradient that is neglected in con- 
ventional hydrodynamics. We will retain this term here 
to make our comparison with the results of the boson 
theory more transparent. Also this term will be needed 
below to provide a large wave vector cutoff to the inte- 
grals determining the renormalized tilt modulus. 



The long wavelength part of the tilt-tilt auto correlator 
can now be evaluated using the definition, Eq. (2.27). 
The non-Gaussian terms in the local hydrodynamic free 
energy ( |6.11 ) are separated out by writing 



F = Fi 



SF e 



3.12) 



where Fq is given by Eq. (4.6), but with the values spec- 
ified in Eqs. (6.1C) for the elastic constants, and 



SF e 



2n Q 



, 8n 



(6.13) 



The tilt autocorrelator is then evaluated perturbatively 
in the non-Gaussian part SF e of the free energy. The 
perturbation expansion is outlined in Appendix B. To 
leading order, we obtain 



^•(q)=TP.(q)+^(q) ! 



(6.14) 



wher e T® ( q) is the bare part of the correlator, given by 
Eq. ( 4.10| - |4.12| ). The hydrodynamic limit of the correc- 
tion STij (q) is given by 



lim £TL(0. 



n k B T , 



[k B Tf 
e\ 2 LA 



E 



(e B ( q < ± )/(k B T))*- q <\ 
[(e B (q> ± )/(k B T))* + q >l 



(6.15) 



This result is identical to that obtained obtained by ular, th e lon g wavelength tilt modulus de fined according 
Tauber and Nelson via the boson formalism. In partic- to Eq. (3.1) is found to be given by Eq. (5.20), with 



B _ n k B T 



E 



{e B ( qi _)/{k B T)f 



2LA ^^[(e B ( q± )/(k B T)r + q 



212 ' 



(6.16) 



which becomes identical to Eq. ( |5.2l| ) in the thermodynamic limit of large sample size. 



VII. TILT MODULUS FROM NONLOCAL, 
NON-GAUSSIAN HYDRODYNAMICS 
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As discussed in the Introduction, neglecting the inter- 
action among vortex segments at different "heights" z 



has severe effects on the flux liquid tilt modulus, namely 
it completely neglects its compressional part, which is the 
largest contribution over a wide part of the (H, T) phase 
diagram. Hence our desire to develop a simple formalism 
for the calculation of the tilt modulus of a flux-line liquid 
that incorporates such nonlocalities. 

A generalization of the boson mapping that incorpo- 
rates the z-nonlocality of the vortex interaction was pea- 
posed some time ago by Feigel'man and collaborators^. 
The z-nonlocality yields a retarded interaction among 
the bosons that can be handled by the introduction of a 
Chern-Simons gauge field. In the limit of infinite pene- 
tration depth, Aj_, considered by these authors, the flux- 
line array then maps onto a charged superfluid. These 
authors argued that the charged boson system possesses 
a normal-fluid phase at zero temperature, corresponding 
to a thermodynamically distinct disentangled flux liquid 

'-sjj(X{[^'-^ + «'-< 



phase, with infinite tilt modulus and longitudinal super- 
conductivity. 

Nonlocality is incorporated in a natural way in hy- 
drodynamics. A nonlinear hydrodynamic free energy 
that incorporates all nonlocalities of the intervortex in- 
teraction can be obtained phenomenologically by coarse- 
graining of the microscopic energy of the vortex liquid, 
following the methods described in Ref. Care must 
be taken in handling the self-interaction between seg- 
ments of the same flux-line at different z heights, which 
is responsible for the non-Gaussian terms in the hydro- 
dynamic free energy. Such non-Gaussian terms are ne- 
glected in the linearized theory, but as seen in the pre- 
vious section they control the renormalization of the tilt 
modulus. The nonlinear hydrodynamic free energy ob- 
tained by such a procedure is given by 

F(r)t H (r') + B{r-r')Sh H (r)Sh H {r') \, (7.1) 



where B{r) is the real space compressional modulus and 
K c {v) is the interaction part o f th e real space tilt mod- 



ulus. The first term in Eq. (7.1) arises from the self- 



energy part of the interaction and it represents a sort 
of nonlinear "kinetic" contribution to the total energy of 
the flux-line array. To make contact with conventional 
notation, it is convenient to rewrite the interaction part 
of the free energy in wave- vector space, 



F 



' f- [t H (r)] 2 



i*(r) 



2nln 



£{c£(q)|F(q)| s 



(7.2) 
c\M)\5n H {ci)\ 2 } 



where the bare compressional modulus, Cj ) 1 (q), and the 
interacti on p art of the bare tilt modulus, cfKq), are given 
in Eqs. ([O]) and (1.3), respectively. 

The non-Gaussian hydrodynamic free energy can also 
be derived from the action of two-dimensional bosons 
with retarded interaction written down by Feigel'man 
and collaborators by successively eliminating nonhydro- 
dynamic fields in favor of hydrodynamic fields via formal 
manipulations analogous to those described in the pre- 
vious section. This derivation is outlined in Appendix 
A. The resulting free energ y di ffers from the phenomeno- 
logical one given in Eq. ( |7.2| ) only in that it contains 
an additional term proportional to density gradients (see 
Appendix A). This term is usually neglected in hydro- 
dynamics because it is of higher order in the gradients. 
We will, however, retain it here as it provides an intrinsic 
large- wave- vector cutoff to the integrals determining the 
renormalized tilt m odu lus. It can be incorporated in the 
free energy of Eq. ( |7.2| ) by the replacement 

c?i(q) -> c?i(q) + (k B T) 2 n q 2 J(4h). (7.3) 



It is convenient for the following to separate out the non- 
Gaussian part of the hydrodynamic free energy of Eq. 
(fF72| ) by letting 



F = F G 



■SF, 



where Fq is given by Eq. (4.6), and 



SF 



I f £i[t H (r)] 2 Sn H (r) 



no 



n H (r) 



(7.4) 



(7.5) 



The tilt autocorrelator can be evalu ated by treating the 
non-Gaussian part of the free energy (7.5) perturbatively. 
Some details are given in Appendix B. The dimension- 
less parameter that controls the expansion in SF/ksT is 
proportional to (ei/2fcsTy / rio) 2 = (/ z / flo) 2 , with l z the 
entanglement length given in Eq. ( |6.6| ). Small values of 
lz/a>o correspond to an entangled flux-line liquid. The 
"kinetic" nonlinearities that are incorporated perturba- 
tively stiffen the tilt modulus of the line liquid, making 
it therefore less entangled. 

The nonlinearities embodied in SF yield corrections to 
all the correlation functions. Here, we only display the 
result for the transverse part of the tilt-tilt correlator, 
that determines the wave vector-dependent tilt modulus. 
Using Eq. (2.30), the wave vector-dependent tilt modu- 
lus is given by 



1 



1 



with 



n ei n n (q±,q z ) 



(7.6) 
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n n (q±,q z ) 



k B T 
LA 



E 



n eik B T 



/2 

9 x 



n ei(q± - q' ± ) 2 



c 4 4 (q') «f* + [6(q')]- 2 c0 4 (q')c° 4 (q- q') (?, - t^) 2 + &(q-q')]- 



E 



(q± • ql) 2 (q± - q'x) 2 fe(q')]" 2 - [l - (4i ■ q'x) 2 ] o'Iq'.W. ^ 



LA *a> C 44(q')^ 4 (q - q') + Ma')]- 2 ] [(<?* - £) a + K.(q - q')]- 2 ] 



(7.7) 



and 



c^(q) 



(fcBT) 2 n g 2 



4ci 



(7.8) 



The length scale £ z (q) differs from the one defined in Eq. (5.19) in that it contains an additional term arising from the 
coupling to the density gradient contained in our free energy and usually neglected in hydrodynamics. For simplicity, 
we use, however, the same notation as in Eq. (5.19). 

The long- wavelength tilt modulus is determined by n„ = lim^^o lim ?i _>o n n (q±, q z ), given by 



k B T 



\ " 



LA ^ c° 4 (q) 



1 - 



i 



n eik B T 



E 



q 2 + [£,(q)]- 2 2LA ^ [ql + Kz (q)]- 2 ] 



(7.9) 



Equations ( |7.6| - |7.9| ) are the central result of this paper. 
If the z-nonloc ality of the intervortex interaction is ne- 
glected in Eq. ( |7.7|) by replacing the elastic constants on 
the right-hand side with the corresponding values u sed i n 
the the lo cal boson formalism, according to Eq. ( 3.10 ), 
then Eq. (7.7) becomes identical to the result obtained by 
TN. In particular, the first term on the right hand side of 
Eq. ( 7.9) is absent in the local boson model of TN, where 
c 44 = n o^i- The long- wavelength normal fluid density is 
then given by Eq. (5.21) and vanishes f or L — > oo. 

The normal fluid density given in Eq. ( |7.9D can be eval- 
uated explicitly for the case of an isotropic superconduc- 
tor (p = 1) in the lim it of infinite thickness (L — ► oo). Af- 
ter inserting in Eq. (7.9) the expressi on f or the nonloca l 
bare elastic constants given in Eqs. (4.7) and ( p~l| - |1.3| ), 
the q z integral in Eq. ([7^) can be evaluated. The re- 
sulting normal-fluid fraction depends on the three length 
scales that characterize the system. These are the aver- 
age intervortex spacing, ao — 1/y/no, the the ab plane 
London penetration depth, A_l , and the single- vortex en- 
tanglement length, l z . We have introduced two dimen- 



sionless parameters, 
2l z 



2ci 



and a dimensionless volume fraction of vortex lines 

1 



(7.10) 



47moAj_ 



(7.11) 



The renormalized long- wavelength tilt modulus is written 
in terms of our dimensionless parameters as 



1 



1 



1 - 



1 + v* no 



and the normal fluid fraction is given by 
n 1 f°° 

— = — dx{K (x\u, v*) + L(x\u, v*)} 
n 2u J 



where 



K (x\u, v* 



x 2 [l + {x + v*)(l + x/u 2 )} + 2z x z 2 x(x + v*) 



yjx + V*Z\Z 2 (z\ + Z 2 )[^/l + X + V*(x + Z\Z 2 ) + Z\Z 2 (z\ + Z 2 ) 



(7.12) 



(7.13) 



(7.14) 



L{x\u, v*) = v 



x(z 1 +z 2 ) 
z\z 2 (z\ + z 2 )(z\ - z 2 ) 



with 



z 1>2 = + x + (x/u) 2 + v*± [(1 + (x/u) 2 -x- v*) 2 + 4v*} 1/2 } 1/2 

v2 



(7.15) 



(7.16) 
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These integrals have been evaluated numerically. The 
resulting normal fluid fraction is shown in Fig. 1 as a 
function of u for several values of the volume fraction 
v*. We note that the dependence on v* is rather weak, 
particularly for small values of it. 




Fig. 1. The normal-fluid fraction given by Eq. (7.13) as a 
function of u for five different values of v* . Notice the 
weak dependence of n n /no on v* for small values of u. 

For v* — (which can be interpreted as either the high 
density limit or the-infinite Aj^ case treated by Feigel'man 
and c ollab oratoral3) , the normal- fluid density given in 
Eq. ( 7.13| ) reduces - up to an overall f^tef of 2 - to 
t he result obtained by Feigel'man et aLEJS Our Eq. 
( (7.9|) generalizes the result obtained by Feigel'man and 
coworkers to the case of finite penetration depth. 

We stress that our calculation is perturbative and we 
have only evaluated the leading correction in perturba- 
tion theory. As discussed above, the small parameter in 
the perturbation theory is proportional to u 2 ~ (h/ao) 2 - 
In other words, the unperturbed state is an entangled 
flux liquid, with a very small value of the z-axis coher- 
ence length l z , and interactions stiffen the vortices, en- 
hancing the tilt modulus. We can estimate the values 
of magnetic field and temperature where our perturba- 
tion theory breaks down as determined by the root of the 
equation 



1 + v* no 



= 1. 



(7.17) 



The solution u (v*) of Eq. ( [7 . 1 7| ) defines a line B D0 (T) 
in the (H, T) phase diagram that can be interpreted as an 
estimate of the phase boundary between entangled and 
disentangled liquid regions. For B > Bdq(T) the liquid is 
entangled, while for B < Bdo(T) the perturbation theory 
breaks down, signaling the appearance of a disentangled 
flux-line liquid. Of course, in order to interpret the region 
B < Br>o(T) as a disentangled flux liquid the Buo{T) line 
must lie in the molten region of the (H, T) phase diagram. 



At high density, v* << 1 and Eq. (7.17) can be approx- 
imated as n n /n ~ 1/v* >> 1. It is clear from Fig. 1 
that the roots of this equation occur at large values of u, 
where n n /no ~ (l/2)ln(u). We then estimate that our 
perturbation theory breaks down for uq(v*) ~ exp(2/u*). 
Converting to field and temperature, this cor resp onds to 
B D o(T) ~ (H c i/2lnK)ln(H c i^o/^k B Tip 2 V^K), with 
H c i = 4>o/A'nX 1 hxK. Below this line, C44 is strongly 
renormalized upward by interactions and a large dis- 
entangled flux-line liquid fraction may appea r. Con- 
versely, at low density, v* » 1 and Eq. ( 7.17 ) becomes 
n n /no ~ 1- The solution of this equation depends weakly 
on v*, as seen from Fig. 1, and is approximately uq ~ 2, 
corresponding to B D0 (T) ~ (0 o /47r)(ei/A: B T) 2 . This re- 
sult coincides with the estimate obtained by Feigel'man 
et al££l, but it applies in a different field regime. The 
solution u Q (v*) of Eq. ( |7.17| ) for general values of v* 
has been obtained numerically and is shown in Fig. 2 as 
a solid line. For small v* (high vortex-line density) Eq. 



(7.17) predicts that the perturbation theory breaks down 
at very large values of u, in a region that is well beyond 
its range of applicability. 



40.0 



30.0 - 



=1 20.0 - 



u,K 



10.0 



0.0 



0.0 



2.0 



4.0 

V* 



6.0 



8.0 



Fig. 2. The solid line is the numerical solution of Eq. 
(7.17). It defines the line uq(v*) in the (u,v*) param- 
eter space where the perturbation expansion of the tilt 
autocorrelator breaks down. The dashed line is ui(v*), 
where n n /riQ = 1. 

We now wish to compare our perturbative result to the 
nonperturbative expression for C44 proposed by Larkin 
and Vinokur and given in Eq. ( |l.6[) . As discussed in the 
Introduction, if the Larkin- Vinokur formula is expanded 
for small values of the normal fluid fracti on n n /no, the 
leading term has the form given in Eq. (1.8), which is 



identical to the long wavelength limi t of our result (|7.q), 
provided we identify n n in Eq. (1.8) with our perturba- 
tive expression for the normal fluid density given in Eq. 
(7.9). It is then tempting to conjecture that a nonpertur- 
bative generalization of our calculation may indeed yield 
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the expression ( |l.6| ) proposed by Larkin and Vinokur for 
the renormalized long wavelength tilt mo dulus, but with 
a normal fluid fraction given by Eq. (7.9), corresponding 
to 



1 



(7.18) 



i/l>0 



with n n given by Eq. (7.£). We stress that Eq. (7.18), 
which is simply a rewriting of the Larkin- Vinokur result, 
is purely a conjecture in the context of our work. It is, 
however, in teres ting to explore its consequences. Accord- 
ing to Eq. ( 7.18 ) , the condition for the vanishing of 1 /c 44 , 
corresponding to the onset of a macroscopic disentangled 
fluid fraction, would read 



n 



1 



(7.19) 



The numerical solution of this equation, denoted by 
ui(v*), is shown in Fig. 2 as a dashed line. We note that 
the line u (v*), where the perturbation theory breaks 
down, and the line Ui(v*), where the conjectured nonper- 
turbative form of l/c^ vanishes, coincide at large v* , but 
diverge at small v* . In this high density region it appears 
that the perturbation theory strongly underestimates the 
stiffening of C44 from interactions. The line ui(v*) defines 
a second "disentanglement line", Br>i(T), in the (H,T) 
phase diagram. Assuming ui(v*) ~ 2 ~ constant over 
the range of v* values of interest, we estimate Box(T) ~ 
(</> /47r)(?i/fc i3 r) 2 . Notice that the field B D1 (T) (which 
coincides with Bdq(T) at low vortex density) is of the 
order of the melting field B m (T) of the vortex lattice. 
Using a Lindeman criterion for melting, this is found to 
be B m (T) = (16ci^ p 2 Mnn) 2 (h/k B T) 2 , where c L is 
the Lindeman parameters. 

Before discussing the location of the disentanglement 
lines Bdq(T) and Boi{T) in the (H,T) phase diagram, 
we recall that the explicit evaluation of the integrals de- 
termining the normal fluid density has been carried out 
for isotropic superconductors (p — 1). To estimate the 
relevance of our result to the anisotropic Cu02 materi- 
als, we have used the above estimate for the boundary 
between disentangled and entangled liquid regions and 
inserted parameter values typical of these materials. To 
justify this approximation, we note that for p >> 1 the 
compressional part of the tilt modulus arising from the 
nonlocality of the vortex interaction in the z direction be- 
comes less important relative to the vortex part. As it is 
precisely this nonlocality that is responsible for a nonva- 
nishing renormalization of C44 in infinitely thick samples, 
we expect that the results that we have obtained for the 
isotropic case will provide an upper bound for the size 
of the renormalization in anisotropic materials. A sketch 
of a phase diagram showing the location of the disentan- 
glement lines Bd (T) (dashed line) and Bui(T) (dotted 
line) is shown in Fig. 3. It is not drawn to scale. 

Using parameter values of YBCO and BSCCO we have 
estimated that in both these materials at high fields 



(B > lTesla) the Bno(T) boundary defining the break- 
ing down of our perturbation theory lies well within the 
flux lattice phase. At low fields there is a possibility for 
a disentangled phase in the reentrant liquid region. This 
region is, however, rather narrow, particularly in YBCO 
wherepit is expected to have a width of the order of 1 
GaussEd. For this reason, while we have drwan in Fig. 
3 the "horizontal" part of the Bno(T) curve as passing 
through this reentrant liquid phase, it could very well 
be that this line is located either above (in the lattice) 
or below (in the Meissner phase) the sketched position. 
The disentanglement line Bdi(T) is shown as dotted in 
Fig. 3 and it is estimated to lie in the liquid phase. The 
existence of this line is, however, just a conjecture in the 
context of our work, as our results are strictly perturba- 
tive. In general we expect the actual disentanglement line 
to lie between our perturbative estimate Boo{T) and the 
conjectured But (T). It could therefore lie almost entirely 
in the solid phase, indicating that a true thermodynamic 
disentangled liquid phase does not exist. This conclu- 
sion would .appear to agree with the latest results from 
simulationaiBcH. Further work beyond the naive lowest 
order perturbation expansion discussed here is needed, 
however, to settle this point. 
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Fig. 3. A sketch (not to scale) of the phase diagram show- 
ing the location of the "disentanglement" lines discussed 
in the text. The dashed line, Bdo(T), marks the break- 
ing down of the perturbation expansion for the inverse 
tilt modulus; the dotted line, Bm{T), corresponds to 
fin/no = 1 and signals the divergence of the conjectured 
form of cf 4 , given in Eq. (7.18). The width of the reen- 
trant liquid phase is in reality much smaller than shown 
here and the line Bdo{T) may or may not pass through 
it. B m (T) is the melting line. H C 2(T) marks the onset 
of a Meissner effect and is not a sharp phase transition. 

One important outcome of our work is that the nonlo- 
cality of the intervortex interaction in the field direction 
has important qualitative effects on the tilt modulus. In 
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particular, it always yields a finite - although often small 
- upward renormalization of C44 even in infinitely thick 
samples. This renormalization is.-absent in calculations 
based on the local boson mappingca. In fact, in the work 
of TN an important role is played by the invariance of 
the flux-line interaction under an affine transformation 
or uniform tilt (corresponding to Galilean invariance of 
a pure boson system). L.D. LandarO has shown that 
the Galilean invariance implies that the superfluid den- 
sity at the ground state (T = 0) of a superfluid equals 
the total density. The affine transformation invariance 



is not present in the more general intervortex free en- 
ergy that allows for pairwise interaction among vortex 
segments at different heights, z. This nonlocality breaks 
the "Galilean invariance" and yields a tilt-tilt interac- 
tion which penalizes any misalignment of the flux-lines, 
therefore favoring disentanglement. 
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VIII. APPENDIX A - DERIVATION OF NONLOCAL HYDRODYNAMICS FROM THE PARTITION 

FUNCTION OF 2D CHARGED BOSONS 



In this appendix we show that the nonlocal, non-Gaussian hydrodynamic free energy given in Eq. (7.2) can be 
derived by fopnal manipulations of the partition function of a two-dimensional charged boson fluid. Feigel'man and 
collaboratorst^l have shown that the partition funct ion of an array of flux-lines described in the London approximation 



by the Ginzburg-Landau free energy of Eq. ( 2.18 ) can be mapped onto that of a two-dimensional system of bosons 
interacting via a massive vector potential. The nonlocality of the intervortex interaction is incorporated via a gauge 
field that mediates a retarded interaction among the bosons. The coherent-state formulation of the boson problem 
yields the imaginary-time action: 

S c [i(j,ip*,a,A] = J drj dr ± { i/>*[ W T + ia - — (fiVi + i&±) 2 - n }ip + V sr {^*) + ^ (Vj_ x a ± ) 2 + 

[z x (d T a ± - Va )] 2 + x a) • A + i-(V X A) 2 } . (8.1) 

The correspondence between vortex and boson variables is summarized in Eq. (|5.4| ). The coupling constant g 
corresponds to the strength of the vortex interaction, according to g 2 <-*■ </)q/(47tA 2 l ) andp is the anisotropy parameter 
that here allows for different scalar an transverse interaction among the bosons. A is the vector potential of the real 
magnetic field (V x A = B), and a = (do, aj_) is a gauge field that mediates the non- instantaneous interaction among 
the bosons. The boson chemical potential \i has to be determined so that the equilibrium boson density Hb equals 
the vortex density, ub = n>o = B/4>o- Finally, V sr is a short range repulsion (on scale £) between the bosons. This 
action is based on the gauge V ■ A = and V_l ■ aj_ = 0. The choice of Vj_ • a^ = instead of V • a = reflects 
the assumption of nonrelativistic velocities for the bosons, corresponding to small tilt of the flux lines aw ay from the 
z directioncll. By rewriting the boson fields in terms of an amplitude and a phase, as defined in Eq. (5.7), we obtain 



fP h f h 2 (V^n) 2 n h 

S c [n, 9, a, A] = / dr dr±\ ihhd T 9 + inao + h V sr (n) + - — a 2 , H n(V0) ■ aj^ — fin 

J Q J 8m n 2m m 



^n(V ± 6) 2 + J* (V.L x a± ) 2 + l[zx (d ra± V« )] 2 



— L—(y x a) • A + -L(V x A) 2 } . (8.2) 
2^/TrX^g Sn 

The assumption of small fluctuations allows us to extend the range of 9 from [— 7r,7r] to [—00, +00]. As described in 
section V, we now eliminate the phase 9 in favor of a vector field P via a Hubbard-Stratonovich transformation, to 
obtain 

S/r- - .1 f 0% , f . r K 2 .aj n 2 . A h 2 (V ± n) 2 Tr ... 

5'n,P,a,A = / dr / dr ± { —nP 2 + — aj + ma Q + — t fj,n + V sr (n) + 

Jo J 2m 2m 2m n 

[l(Vixa 1 ) 2 + i[zx(9 r a 1 -Vflo)] 2 + n -L (V x a) ■ A + 

1 , . ,9 nnh 2 , , n . 1 . . 

Vx A) 2 + ^ln( — ) }, (8.3) 
87r m no 
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71 

d T n + V ± -—(hP + isi ± ) = 0. (8.4) 
m 



with the constraint 



The last term in the action in Eq. ( p. 3D , logarithmic in the density, is the Jacobian of the transformation. We then 
make a change of variables, 

a 77, * 

t=— (ftP+iax), (8.5) 
m 

and obtain 

S c [ra,t,a,AJ = / rfr / dr ± { — ?a_i_ • t + ma - fm + V sr (n) + h 

Jq J Zifl oTYl Ti 

^ (Vi x aj 2 + l[z x (3 T a ± - Va )] 2 + -^(V x a) ■ A + i(V x A) 2 }, (8.6) 

with the constraint 

d T n + V_l • t = . (8.7) 

The Jacobian of this transformation cancels that of the previous one. 

Finally, we define an effective action S° s for the bosons by integrating out both the vector potential A(r) and the 
gauge field a(r), 

J VnVtVAVa e-^'^'^^n + V_l • t) = J VnVte- s ^ [h ^8{d T n + V_l ■ t). (8.8) 

The prime over the integral sign on the left hand side of the equation indicates that the integration over A and a 
has to be performed by taking into account the constraints imposed by our choice of gaug e. The vector potential and 



gauge field are most easily integrated out by rewriting the field part of the action (3.6) in Fourier space, with the 
result, 

Sf[n, t] = dr [dr ± {^-^n + V sr (h) + -^^^ } 

+ —Y{ _lt T (q)] 2 + 4 g2X \ In(q)l 2 }, (8.9) 



where ty(q) = x t(q). By making use of the continuity constraint given in Eq. (8.7), we can write 



n 2 \ 2 n 2 n 2 \ 2 n 2 \ 2 

9 K |iT(q)| 2 + \-^4-Rq)| 2 = : 2 \ ^ It(q)| 2 



1 + q 2 X\ + q]_ p 2 X 2 L Q±l + q 2 X ± 1 + £?f + q\ p 2 X\ 

g 2 X 2 (l + q 2 p 2 X 2 L ) 



|n(q)| 2 . (8.10) 



(l + q 2 Xl)(l + q 2 Xl + qlp 2 XlY 
Finally, if we replace the short range repulsion V sr (n) by a short-wavelength cutoff and identify the boson density h 



and momentum field t with the corresponding hydrodynamic quantities for the vortices, we see that Eq. (8.9) yields 
precisely the nonlocal non-Gaussian hydrodynamic free energy discussed in section VI. 

IX. APPENDIX B - PERTURB ATIVE CORRECTIONS TO THE TILT MODULUS FROM NONLINEAR 

HYDRODYNAMICS 

The wave-vector dependent tilt modulus is defined in terms of the transverse part of the tilt-tilt correlator as in 



Eq. (3.1). In the hydrodynamic approximation, the tilt-tilt correlator can be written as 

jVh{Y)Vi{v)k{Y)ij{v')e- F / kBT 8(d z h + Vj_ • t) 
/M(r)Dt(r)e- F / k BT^a z n + Vj_ • t) 
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Where F is the hydrodynamic free energy given in Eq. (|7.2|). The free energy can be written as the sum of a Gaussian 
part, Fq, and non-Gaussian corrections, 6F, as in Eq. fl7.4| ). We want to calculate up to lowest-order in the small 
parameter, it 2 , nonlinear corrections to the tilt autocorrelator. By keeping only terms up to fourth order in the 
fluctuations of the hydrodynamic fields, the non-Gaussian part of the free energy is given by, 



SF: 



2n 2 ft 2 



*i( t li)*i( t l2)^n(-qi -q 2 ) 



qi,q2 



2ngft 3 



Y *i( t li)*i( t l2)^(q3)^n(-qi - q 2 q 3 ) • 



(9.2) 



qi ,q2 ,q3 



The tilt-tilt correlator is then evaluated in Fourier space perturbatively in the non-Gaussian part of the free energy, 
with the result, 



'G > 



(9.3) 



where (•••)£; denotes a cumula nt a verage over the Gaussian ensemble wi th weight ~ cxp(— Fa/hsT). The first term 
on the right hand side of Eq. ( |9.3| ) is the Gaussian result given in Eqs. (4.10-4.12). 

Using Wick's theorem, the corrections arising from the non-Ga ussian pa rt of the free energy are easily expressed in 
terms of the correlations in the Gaussian ensemble given in Eq. (4.8-4.12), with the result 
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and 



P^(qx)^(q)^(q')(^) 2 )& = 2n<Wo[T£(q)] 2 i£(q±)^ ]T {(if (qi)if (- qi )) G (\Sn H (q - qi )| 2 ) G + 

qi 

(if ( qi )^(- qi )) G (f (qj - cOSh H (q- qi )) G }. (9.5) 

By substituting the expressions for the Gaussian correlators given in Eqs. (4~§|-4.12), we obtain the following expression 
for the transverse part of the tilt autocorrelator to lowest order in the non-Gaussian terms, 
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